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Abstract 



We consider the inclusive cross section for jet production with large trans- 
verse momentum in deep-inelastic scattering. This process has been proposed 
as a probe of small-x physics, particularly the measurement of 'hot spots' in- 
side the proton. We present a numerical calculation of this process, taking 
into account a larger phase space. The theoretical reliability as well as phe- 
nomenological uncertainties of the calculation are discussed. 
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1 Introduction. 



In order to study the main properties of low xb deep- inelastic scattering 
processes, Mueller suggested in ^ an experiment in which all anticipated 
new phenomena in this kinematical region should have a large effect. His 
idea was to measure the inclusive production of a gluon jet with a transverse 
momentum kj^t very close to the photon virtuality Q and with a fraction of 
energy Xj as close to one as is feasible, so that the ratio xb/xj can be small. 
In this case, 

• the cross section of the process can be calculated within the framework 
of perturbative QCD, if ^ >> Aq^^.^. 

• the dependence of the cross section on is governed by low xb gluon 
emission which can be described by the BFKL evolution equation 
This is in contrast with the usual GLAP approach, in which the 
cross section is described by the simple Born diagram of Fig. la. and 
turns out to be constant with respect to x^. 

• the scale of the shadowing corrections is determined by the size of the 
'hot spot' , namely R ~ l/kjt and they are expected to be large (see 
ref. i). 

Numerical estimates for this process have been performed in a series of 
papers 0-0, but, in our view, the matter has not been settled yet (cf. the 
strong dependence on an infrared cut-off in the BFKL equation with running 
coupling, for values of k"^ and that are smaller than about 50 GeV^ 

In this paper we reconsider the theoretical formulae for the cross section 
and present our numerical estimates in studying the small-x and infrared 
behavior of the one-jet inclusive cross section. We will not consider any 
shadowing corrections. 



2 The basic formula. 

The correct formula for inclusive one-jet production in the region of small xb 



was given in ref. [|TT[. In terms of a differential structure function involving 
the jet variables xj and kf (resp. the longitudinal momentum fraction of the 
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jet and its tranverse momentum squared), it looks as follows for gluon jet 
production 



- ■ asi'min{k.t,kot,kt}) 



b{—, kl, Q')Mxj, kl,) ■ 5(2) (kt - ki,t - k2,t) 



3Nc f dkl^dkl^ ■ r,2 ,2 ;2in 

^kV ^-Xikl, kl, kj) 

kl,Q')<Pc{x,,kl), (2.1) 



Xj 



where all notation is explained in Fig. lb., X{x, y, z) = x'^ + y"^ + — 2xy — 
2xz — 2yz, and (pc is the gluon density, related to the gluon distribution 
function xG{x, Q"^) by 

das(k'^)xG(x,k'^) , , ,0, ,,0, ^. 
^ ^ ' = k^)as{k^). (2.2) 

As was shown in ref. |jlO|, eq. (|2.1| ) can be reduced in the double logarithm 
approximation (DLA) of perturbative QCD to the expression of eq.(51) in 



n|], where the inclusive production in deep-inelastic scattering was studied 



in detail within this approximation, (see also [|l^). If we integrate in eq. 
( |2.1| ) only the part of phase space where k^^ « k\^ [kf^ ^ kf) ('small' 
phase space) we can rewrite eq. (|2.1j ) in the form: 



dF2 3N^ 



d In Xidkf 71 k? 



/ dkl t ■ as{kl t)(pB{ — , k^, Q'^)(f)G{xj, kl ^) 

J Xj 

asik'M — , kl g2)x,G(x,, k^. (2.3) 



It is this equation that was used in all previous numerical estimates of the 
inclusive jet production Here we calculate the differential structure 

function without any restriction on the region of integration ('large' phase 
space). In comparing the two we will find an enhancement of the cross section 
due to the larger phase space of about 80 %. 
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3 Calculational procedure. 



The functions in eq. (2J.) are solutions of the BFKL-equation 0, 



ain(l/x) ~ ~ir kl ^ \k'^-k^\ F2(4fc'4 + fc4)(i/2)/- 

(3.1) 

We now discuss aspects of our procedure of solving this equation. 

We chose different initial conditions for the functions 0_b(^, k'f,Q'^) and 

(pci^j, kl) in (|2.1|). For 0_b(^, k^, Q"^) we used the same initial condition as 

in ref. 0, namely aX zq = xs/xj = 10^^ 

2 2^ F,{zo,kj,Q^) Fo{e,Q') 
0b(2;o, K,q) = p ~ p , (3.2) 

where the function Fq, related to the quark box diagram, was calculated in 
refs. For (pd^j, k^) we have to reconstruct the initial condition from 

experimental data. However, in order to solve the BFKL-equation we need 
to know the behavior of cf) at any value of k'j for some fixed x, even at k'j 0. 
To accomplish this, we used the following procedure to describe the low k^ 
behavior of XjG{xj, kl): 

2 k"^ ^ ^ ~ 

k^ + 



where k'^ = k"^ + q^, x = x/{x + k"^ /k'^{l — x)). We used the mapping (|3T 
because (1) this parametrization ensures XjG{xj, k'^) ~ /c^ as A;^ — and 
such a behavior is the direct consequence of the gauge invariance of QCD, 



and (2) it works for the case of F2 (see refs. |T^-[| 



For the function XjG{xj, k"^) we used a fit Q to the data set from the CTEQ 



collaboration []14| down to Q values of 1 GeV . The initial condition for (j)Q 
was then constructed according to ( p.2| ) (we used a fixed coupling), at a value 
X = 10~^, where, as in 0, instead of xG{x, fc^) we used the effective density 
xG{x, k^) + I E/=i x[qf{x, k^) + qf{x, k% 

We solved^ the BFKL equation ([Ol ) with a fixed coupling 0:5 = as{Q^) 
because this equation only sums the leading logs (0:5 ln(l/a;))", and not the 

^We are grateful to J. Botts for making this fit, and discussions relating to it. 
^We are grateful to J. Kwiecinski for sending us his program to solve the BFKL equation 
numerically. 
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subleading ones. The latter remains an unsolved problem. However, we made 
a rough estimate of how important the corrections from a running as could 
be by calculating an average as 

The dependence of < as >i on xb showed that we could use as = CisiQ'^) as 
a good first approximation. For example, for k"^ = Q"^ = 10 GeV^, we found 
asiQ'^) ~ 0.2, and < as >i~ 0.15 for x < 10"^ for both i = B,G. 
A much more detailed study of this problem can be found in 
The last point we discuss is the infrared cut-off in eq. ( p.l| ) . In principle 



eq. ( |3.1|) is infrared stable and one can take /cg = 0. However, in order to 
further investigate the dependence of the solution of ( |3.1D on small momenta 
we preferred to introduce a cut-off k^ and see how much the answer depends 
on its value. 



4 Results and conclusions. 

The results are shown in Figs. 2a and 2b. for the values k"^ = Q"^ = 10 GeV^, 
Xj = 10~^. The first observation we make from Fig. 2a is that the answer for 
( |2.1|) when one includes the correct ('large') phase space increases the results 
of previous calculations (small phase space) by about 80 %. Furthermore, 
the dependence on is visible, but does not compensate for the difference 
between small and large phase space. A similar conclusion we found to hold 
for the dependence on the infrared cut-off k^. 

The most discouraging result is shown in Fig. 2b. Here we plot the nor- 
malized differential structure function 

o 1 dF2{xB,Q^;Xj,k^) 

^ F2{xb,Q^) dlnxjdk^ ^ ' ' 

where we took F2{x,Q'^) = J dk'^cpB^^, k"^ ,Q'^)4>G{xj, k"^). One notes that 
this ratio seems to be independent of x^. This would seem to indicate that, 
within the approximations made, this special environment does not seem to 
be much better for measuring small-x effects than a direct measurement of 
F2. This feature persisted when we took F2{x, Q"^) constructed from the MRS 
D_' distributions |I5| . 



5 



References 

[1] A.H. Mueller, Nucl.Phys.B (Proc. Suppl.) 18C, (1991) 125. 

[2] E.A. Kuraev, L.N. Lipatov, V.S. Fadin, Sov. Phys. JETP 45, (1977) 
199; Ya.Ya.Balitsky, L.N. Lipatov, Sov. J. Nucl. Phys. 28, (1978) 822. 

[3] V.N. Gribov, L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; L.N. 
Lipatov, Yad.Fiz. 20 (1974) 181; Yu.L. Dokshitzer, Sov.Phys. JETP 46 
(1977) 641; G. AltareUi, G. Parisi, Nucl. Phys. B126 (1977) 298. 

[4] L.V. Gribov, E.M. Levin, M.G. Ryskin, Phys. Rep. ClOO (1983) 1. 

[5] J. Bartels, A. De Roeck, M. Loewe, Zeit.Phys. C54, (1992) 635. 

[6] W.-K. Tang, Phys. Lett. B278 (1992) 363. 

[7] J. Kwiecinski, A.D. Martin, P.J. Sutton, Phys. Lett. B287 (1992) 254; 
Phys. Rev. D46 (1992) 971. 

[8] A.J. Askew, J. Kwiecinski, A.D. Martin, P.J. Sutton, Durham preprint 
DTP/92/78, December (1992). 

[9] J. Bartels, H. Lotter, DESY preprint, DESY-93-027, March 1993 

[10] K. Charchula, E. Levin, in 'Physics at HERA', Proceedings of the Work- 
shop, DESY, Hamburg, eds. W. BiichmuUer and G. Ingelman, vol. 2, 
(1991), 223. 

[11] L.V. Gribov, Yu.L. Dokshitzer, S.I. Troyan, V.A. Khoze, Sov. J. Nucl. 
Phys. 88, (1988) 1303. 

[12] B. Badelek, J. Kwiecinski, Phys. Lett. B295, (1992) 263. 

[13] H. Abramowicz, E.M. Levin, A. Levy U. Maor, Phys. Lett. B269, 
(1991), 465. 

[14] J. Botts, J.G. Morfm, J.F. Owens, J. Qiu, W.-K. Tung, H. Weerts, Phys. 
Lett. B304, (1993) 159. 

[15] A.D. Martin, R.G. Roberts, W.J. Stirling, RAL preprint RAL-92-078, 
1992. 



6 



Figure Captions 

Fig. la. Born diagram for single gluon jet production in deep-inelastic scat- 
tering. 

Fig. lb. Single-jet production in deep- inelastic scattering in the large phase 
space case. 

Fig.2a. Plot of ^i^^%k2 vs. xb- The value of Xj is 0.01 and we took k'^ — 
— 10 GeV^. The sohd hne corresponds to the small phase space case, 
the three remaining to the large phase space case with three different 
values of c^, namely Q'o = 1 GeV^ (long-dashed), go — 2 GeV^ (short- 
dashed) and go = 4 GeV^ (dotted). 

Fig.2b. The ratio R2 (4.1) for three values of in the large phase space 
case. The notation is the same as in Fig. 2a. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305341vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305341vl 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305341vl 



This figure "fig4-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9305341vl 



